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Using an exact diagonalization treatment of Ising and Heisenberg model Hamiltonians, we study 
field-induced phase transition for two-dimensional antiferromagnets. For the system of Ising anti- 
ferromagnet the predicted field-induced phase transition is of first order, while for the system of 
Heisenberg antiferromagnet it is the second-order transition. We find from the exact diagonaliza- 
tion calculations that the second-order phase transition (metamagnetism) occurs through a spin-flop 
process as an intermediate step. 

PACS numbers: 75.40. Cx, 75.40.Mg, 75.10.-b, 75.30.-m 
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I. INTRODUCTION 



Much earlier (in, 1932) the antiferromagnetic order was 
proposed by Neelu in order to explain the low tempera- 
ture behavior of the magnetic susceptibility of certain 
metals. Currently physical properties of low-dimensional 
quantum antiferromagnets at low temperature are ac-i 
tively pursued. The exact solution via Bethe AnsatzEl 
is limited to one-dimensional integral systems, but not 
extendible to multidimensional systems. Recently, one 
among the interesting subjects is a numerical study of 
field-induced phase transition, namely, the antiferromag- 
netic to ferromagnetic transition in the two-dimensional 
systems of antiferromagnetically correlated electrons un- 
der externally applied magnetic field. By applying the 
dynamical mean field theory (DMFT) to the Hubbard 
model, Held et alB studied the microscopic origin of 
metamagnetism along an easy axis in antiferromagnets 
under external magnetic field. Their study showed that 
at half filling a metamagnetic phase transition arises 
via the first-order phase transition at low temperature 
and that the second-order phase transition occurs near 
the Neel temperature. On the other hand, Bagehorn 
and HetzeO observed the second-order phase transition 
at zero temperature from the projector quantum Monte 
Carlo (PQMC) calculation of the Hubbard model with 
an easy axis. Earlier using a Landau theory of free en- 
ergy Moriya and Usamit] revealed the second-order phase 
transition at low temperature involving the mixed phase 
of ferro- and antiferromagnetic states. However, Bage- 
horn and Hetzel questioned whether the mixed phase 
would survive at a small region of magnetic field if exact 
electron correlations were implemented. In the present 
study, by performing the exact diagonalization calcula- 
tions of the two-dimensional systems of antiferromagnet- 
ically correlated electrons under applied magnetic field, 
we examine the nature of field-induced phase transition 
(metamagnetism) at K for both Ising and Heisenberg 
antiferromagnets. 



II. LANCZOS EXACT DIAGONALIZATION 
CALCULATIONS OF FIELD-INDUCED PHASE 
TRANSITIONS 

With the inclusion of Zccman coupling term, the two- 
dimensional t-J model Hamiltonian is written, 



H 



+J ^(S, • Sj - ^n.rij) - h^<jc\ a c la , (1) 

(ij) ia 

where Ci a (cj CT ) is the electron annihilation (creation) op- 
erator, Si is the electron spin operator at site i, and h 
is the Zeeman energy, (J, B B. Here \x B is the Bohr mag- 
neton and B is the strength of external magnetic field. 
a = +1 (—1) for up spin (down spin). In the case of half 
filling, the t-J model Hamiltonian effectively reduces to 
the Heisenberg model Hamiltonian. 

The above Hamiltonian in Eq. ([[]) will be diagonal- 
ized by the Lanczos exact diagonalization methodfl For 
the exact diagonalization treatment of antiferromagnetic 
correlations, uniform and staggered magnetizations are 



defined by,l 



<K) 2 > 



(2) 



where t 



Here (w£) 2 represents the square of 



x,y,z. — y.» g 
the uniform and the staggered magnetizations in the £th 
direction corresponding to q — (0, 0) and q — Q = (tt, tt) 
respectively. N is the total number of lattice sites. We 
calculate the variation of the uniform and the staggered 
magnetizations with vertically applied magnetic field to 
a 4 x 4 square lattice of antiferromagnetically correlated 
electrons by satisfying periodic boundary conditions. For 
the parameters we choose t = 1 and J = 0.4. Lanczos 
iteration is terminated when the ground state energy is 
converged within the error bound of 10 -10 . 



1 



For the Ising systems, H — JJ2{ij) SfSji we nn d that 
the metamagnetic phase transition is of first order, as 
is shown in Fig. |} In our previous workjj we calcu- 
lated the staggered and the uniform magnetizations for 
the Ising system by using Hubbard model Hamiltonian 
in mean-field approximation. The first-order phase tran- 
sition was also predicted. In both approaches we note 
that there exists no mixed phase (of both m Q ^ and 
m 7^ 0). In other words, discontinuity from staggered to 
uniform magnetization is observed at a critical magnetic 
field, say, h = 0.4i in Fig. || (the staggered magnetization 
beyond the critical field and the uniform magnetization 
below the same point did not completely vanish owing 
to finite-size effect). Thus in the Ising systems metam- 
agnetism occurs as a first-order phase transition. This, 
is in agreement with the DMFT study of Held el alia 
On the other hand, Bagehorn and Hetzel showed from 
their PQMC calculation that metamagnetism occurs via 
a second-order phase transition. 
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FIG. 1. Staggered and uniform magnetizations for the 
two-dimensional Ising system as a function of Zeeman energy, 
(a) staggered magnetization m Q and (b) uniform magnetiza- 
tion m . 



For the Heisenberg systems, we show the predicted 
magnetization as a function of applied field in Fig. |^. 
The stepwise curves are inevitably formed owing to the 
size effect. As the number of lattice sites increases, the 
stepwise curves are expected to turn smooth. The solid 
curves in the figure are interpolated ones only for the 
visual aid. The second-order field-induced transition is 
observed for metamagnetism at K. The magnetic sus- 
ceptibility, x ( = 7pf) changes from a negative to a positive 
value at a crossing point between the lines of the stag- 
gered magnetization (solid line) and of the uniform mag- 
netization (dotted line). This indicates the occurrence of 
a mixed phase of the fcrro- and antifcrromagnctic phases 
during the process of metamagnetism. The predicted an- 
tiferromagnetic susceptibility vanishes beyond a 'critical 
field' corresponding to the Zeeman energy, h ~ 0.8/ (the 
staggered magnetization did not completely vanish ow- 
ing to the size effect). The uniform magnetization is pre- 
dicted to remain constant beyond the same point. Al- 
though not to be directly compared, due to differences in 
the space dimension of antiferromagnet, such a tendency 
of uniform magnetization was experimentally observed 
for the three-dimensional crystal of Y(Coi_ x Al x )2 at low 
temperature by Goto et aln 
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FIG. 2. Staggered m Q (solid line) and uniform m (dot- 
ted line) magnetizations for the two-dimensional Heisenberg 
system as a function of Zeeman energy. Solid curves are guide 
lines for the eyes. 

In order to find the cause of the second-order phase 
transition with the Heisenberg antiferromagnet, in Fig. || 
we show both the staggered and uniform magnetizations 
in the direction parallel to and perpendicular to the plane 
respectively. As the external magnetic field increases, the 
x- and y-components of the staggered magnetization are 
predicted to increase particularly in the region of low 
field, h ~ 0.3/ in Fig. ||(a), whereas the z-component de- 
creases in the same region. The antiferromagnetic order 
projected onto the plane (perpendicular to the external 
magnetic field) tends to persist to a point, while its z- 
component m z Q monotonically decreases and vanishes at 
h ~ 0.8/:. This can be further explained as follows. In 
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the absence of magnetic field, spins align into a strong 
antiferromagnetic state due to the Heisenberg interac- 
tion as is shown in Fig. [|(a) • As the magnetic field in the 
z-direction increases, spin-flop process occurs as an inter- 
mediate state by exhibiting the components of antiferro- 
magnetic spin alignment on the x-y plane, as depicted in 
Fig. |](b). For the Heisenberg antiferromagnet, the spin 
flop, that is, the x- and y-components of the staggered 
magnetization peak at a particular value of the applied 
magnetic field, say, h ~ 0.3i as is shown in Fig. ||(a). As 
the external magnetic field further increases, the Zecman 
effect begins to dominate the Heisenberg interaction with 
disappearance of the spin-flop process by allowing ferro- 
magnetic configurations, as is shown in Fig. f|(c). On 
the other hand, for the fsing system the first-order phase 
transition was predicted with no involvement of spin-flop 
process before the metamagnetic transition, as is shown 
m Fig. |. Thus we claim that antiferromagnetic spin 
interactions involving x- and y-components are responsi- 
ble for the second-order phase transition accompanying 
a spin-flop process as an intermediate step. 
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FIG. 4. Spin flopping and metamagnetism under ap- 
plied magnetic field. The magnetic field is applied in 
the 2-direction. The horizontal solid line represents the 
two-dimensional x-y plane, (a) antiferromagnetic spin order, 
(b) spin flop, and (c) ferromagnetic spin order. 



III. CONCLUSION 



By applying the exact diagonalization method to 
Heisenberg antiferromagnets under external magnetic 
field, we examined how the second-order phase transi- 
tion occurs from an antiferromagnetic to a ferromagnetic 
state at a critical magnetic field. The second-order phase 
transition occurred in the presence of the mixed phase 
near the critical magnetic field. According to the present 
exact diagonalization study, the second-order phase tran- 
sition was induced through the presence of spin-flop pro- 
cess owing to the influence of the x- and y-components 
of antiferromagnetic spin interactions. In short we ob- 
served that for the case of two-dimensional Ising systems 
the field-induced transition from the antiferromagnetic 
state to the ferromagnetic state is of first order with no 
involvement of spin-flop process, while for the case of 
two-dimensional Heisenberg antiferromagnetic systems it 
is of second order with the involvement of the spin-flop 
process as an intermediate process. 
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FIG. 3. Staggered and uniform magnetizations for the 
two-dimensional Heisenberg system as a function of Zeeman 
energy, (a) Decomposition of staggered magnetization (m Q ) 
and (b) decomposition of uniform magnetization (m ) into 
directions parallel (m^ and nig) and perpendicular (m^) to 
the plane respectively for q — Q, 0. 
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